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1. Introduction 



We prove the global existence of H 2 solutions of the equations of barotropic, compressible magneto- 
^ hydrodynamics (MHD) in three space dimensions: 

pt + div(pu) = 0, (1.1) 
(pu j ) t + div(pu^u) + P{ P ) Xj + (||B| 2 ) X . - div{B j B) = pAu j + XdWu Xj , (1.2) 
B\ '+ div(B j u - u j B) = 0, (1.3) 



H. divB = 0. (1.4) 

& . 

Here x € M 3 is the spatial coordinate, t > is time, p, u — (w 1 ,?/ 2 ,?/ 3 ) and B = (B 1 , B 2 , B 3 ) are 
the unknown functions of x and t representing the density, velocity, and (scaled) magnetic field in a 
compressible ionized fluid, P — P{p) is the pressure, p, and A are positive viscosity constants, div and 
A are the usual spatial divergence and Laplace operators. 

The above system (|l.ip - (|1.4l) is obtained by combining the Navier-Stokes equations for compress- 
ible barotropic flow with Maxwell's equations in free space and the ideal Ohm's law. We briefly outline 
the derivation (see Cabannes [2] for details) Using Ampere's law, we have 

dE 

£o~Qj- + J = Po 1 curlB, 

where E, J, and B are the electric field, current density and magnetic field, and £o,Mo are constants. 
The first term on the left is taken to be negligible in the given application, and so by the Biot-Savart 
law, the magnetic force per volume of fluid is 

, , div( BB T ) (\B\ 2 

J x B = p^ 1 evil B x B = ^ '- - V ^— ^ 

Po \ 2po 
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and it gives the forcing terms following the pressure in (|1.2p after scaling out the factor p,Q, On the 
other hand, we combine the ideal Ohm's law E + u x B = with Faraday's law so that 

dB 

— — = — curLE = curlfu x B), 

and together with the Gauss's law for magnetism (|1.4[) . we obtain equation ()1.3|) . 
The system ()l.lj) - ()1.4|) is solved subject to initial conditions 

(p(;0),u(;0),B(;0)) = { Po ,U ,B ), (1.5) 

where po is bounded above and below away from zero, div Bq = in a suitable sense and modulo 
constants, (po,uq, Bo) is small in i/ 2 (R 3 ). 

The subject of MHD was first initiated by Alfven [1] in 1940's, since then it has become one of 
the most challenging topics in fluid dynamics. In the fully viscous case, we further have the resistivity 



term vB^ on the right side of ()1.3[) . so that the system becomes 

p t + div(pu) = 0, (1.6) 

(pu j ) t + div(pu^u) + P(p) Xj + (±|.B|% - div(^B) = nAu* + Adivuay, (1.7) 

B{ + div{B j u - u j B) = vB\ (1.8) 

di Vj B = 0, (1.9) 



where v > is the resistivity constant. A large variety of solutions to (|1.6p - (|1.9l) can be obtained for 
different initial data stated in (|1.5|) . When (po,uo, Bo) is taken to be close to a constant in ii 3 (R 3 ), 
Kawashima [14] proved global existence of "small-smooth" solutions to the system (|l.l|l - (|1.4|) . and 
such solutions remain close to the same constant in £T 3 (R 3 ). On the other hand, for initial data with 
arbitrary large energy and nonnegative density, Sart [20] and also Hu and Wang [HJ-P3] showed 
global existence of "large-weak" solutions to (ll.l[) - (|1.4p based on a method introduced by Lions [TB] 
and Feireisl [I]- [5]. Solutions in this general class possess very little regularity, and some of those 
solutions may be even non-physical (see [9] and [IDJ). In between those two types of solutions just 
mentioned, there is another type called "intermediate-class" solutions which is first introduced by 
Hoff [6] for Navier-Stokes equations and later extended by Suen and Hoff [23] for MHD. In this case, 
initial data is taken to be small in L 2 (R 3 ), and initial densities are assumed to be nonnegative and 
essentially bounded. Solutions posses just enough regularity which makes it possible for one to develop 
uniqueness and continuous dependence theory (see Hoff and Santos pj] for Navier-Stokes equations). 

On the other hand, for the case when resistivity v becomes zero (system (|l.l|) - ([1.4jl ). there is not 
much known result in the literature. In Kawashima |14j . assuming initial data in Sobolev spaces H s 
with higher regularity (s > 5/2), both local and global existence theory of smooth solutions for the 
Cauchy problem of (|l.ip - (|1.4j) were established. The goal of the present paper is thus to generalize 
Kawashima's result by establishing global existence of smooth solutions with small initial if 2 data. It 
seems to be less trivial than one could expect due to the absence of the damping mechanism in the 
transport equation of B. We use a kind of standard energy estimate similar to that which is used for 
the incompressible viscoelastic fluids suggested by Lei-Liu-Zhou [15]. It enables us to extract some 
weak dissipation from the system which is crucial in our analysis. 

We introduce two variables associated with the system (jl.ip - (jl.4l) which are important to our 
analysis. The first one is the usual vorticity matrix uj = u^ ,k = u J Xk — u Xj , while the other one is the 
effective viscous flux F given by 

F = (p + A)div u - P(p) + P(p). (1.10) 

where p is a positive constant density. By adding and subtracting terms, we can rewrite the momentum 
equation (|1.2p in terms of F and w: 

pu> = F X] + - {\\B\ 2 ) X] + div(B'B). (1.11) 
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The decomposition also implies that 



AF = div( 5 ), 



(1.12) 



where g 3 = pu 3 + (^\B\ 2 ) Xj — div(_B J \B). We refer to Hoff [6] for a more detailed discussion of F. 

We now give a precise formulation of our results. First concerning the pressure P we assume that 

(1.13) P G C([0,oc)) n C 3 ((0,oo)),P(0) = and P(p) > for p > 0; and there are positive densities 
p" > p' such that P is increasing on [0,p'} and on [p",oo), P(p) > P(p') for p > p', and 
P(p) > P{p") for p > p". 

Next we fix a positive reference density p and then choose positive bounding densities p and p satisfying 

p < min{/5, p'} and max{/5, p"} < p, (1-14) 

and finally we define a positive number 6 by 

5 = min{min{p, p'} -p,p- max{p, p"}, |(p - p)}. (1.15) 

(Notice that 5 need not be "small" in the usual sense.) Concerning the diffusion coefficients p, and A, 
we assume that 

p,X>0 (1.16) 
Concerning the initial data (po,ito,-Bo) we assume that there is a positive number d < 8 such 



that 

We assume also that 
and we write 



p + d < ess inf po < ess sup po < p — d. 
divS = 

Co = ||P0 - P\\EP + ||«o||fl> + ||-Bo||fla. 



(1.17) 
(1.18) 
(1.19) 



Weak solutions are defined in the usual way; we say that (p, u, B) is a weak solution of (jl.lj) - 
([L5f provided that (p - p, pu,B) G C*([0, oo); i?" 1 ^ 3 )) with (p,u,5)| t = = (p ,u ,5 ), Vit, V5 G 
L 2 (K 3 x (0,oo)), div5(-,i) = in £>'(R 3 ) for t > 0, and the following identities hold for times 
t% > £i > and C 1 test functions ip having uniformly bounded support in x for ( G [t±, £2]: 

*2 /.t 2 



p(x, -)f{x, -)dx 



(pipt + pu ■ Vtp)dxdt, 



(1.20) 



and 



(pu 3 )(x,-)cp(x,-)dx 



j 

J®? 



B J [x, -)(p(x, -)dx 



to 



[pu J ipt + pv?u ■ V<p + P(p)ip Xj ]dxdt 

[\\B\ 2 ip X] ~ B 3 B ■ V<p] dxdt 
[/iVm j • V<p + A(divu)<p x Jc&z;e££, 

/ / [{B j u - u 3 B) ■ \7ip]dxdt. 



(1.21) 



(1.22) 



We denote the material derivative of a given function v by v = Vt + Vu • tt, and if X is a Banach 
space we will abbreviate X 3 by X when convenient. Finally if / C [0, oo) is an interval, C 1 (I; X) will 
be the elements v € C(I;X) such that the distribution derivative v t G £>'(R 3 x intl) is realized as an 
element of C(P,X). 

The following is the main result of this paper: 
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Theorem 1.1 Assume that the system parameters P,p,\ in (|1.1[) - (|1.4I) satisfy the conditions in 
(1.13V ([1.16P and let 6 be as defined in (|1.15[) . Let positive numbers N and d < S be given. Then 
there are positive constants a, C , and 9 depending on the parameters and assumptions in (1 . 13)- (|1 . 16|) . 
on N, and on a positive lower bound for d, such that, if initial data (po,Uo, Bq) is given satisfying 
(fi~T7)) - (fi~T9l) with 

Co < a, (1.23) 
then there is a solution (p,u,B) to (jl.l|) - ()1.5|) in the sense of (|1 .20(1 .22|) on all o/R 3 x [0, oo). The 
solution satisfies the following: 

d 3 t D%u g L°°([0, oo);H 2 - 2 i-^) n £ 2 ([0, oo); ffS-^-M+i), (1.24) 

dlD« Pl dlD«B g L°°([0,oc);ff 2 - 2 ^H), (1.25) 
/or o/Z J, a satisfying 2j + |a| < 2, tujf/j 

P < /K^, t) < p a.e. on R 3 x [0, oo), (1.26) 

and 

sup (||(p-p,u,S)(-,s)||| f2 +||(p t ,u t ,B f (-,s)||| 2 ) , 

0<s<t 

+ T (IIVh(., s )||^ 2 + IKMH 2 ^) d s < ccl (1.27) 

Jo 



The proof is given below in sections 2-5 and begins with a number of a priori bounds for local- 
in-time smooth solutions. The existence of such smooth solutions is therefore crucial, and for this we 
rely on the following result of Kawashima [14], pp. 34-35 and pp. 52-53: 

Theorem (Kawashima) Assume that p and A are strictly positive and that the pressure P g C 3 ((0, oo)). 
Then given p > p > and C3 > 0, there is a positive time T depending on p, p, and C3 and on the 
system parameters p, A and P, such that, if initial data (po — p, uq, Bq) is given satisfying 

IK/00 - P,U ,Bq)\\ H 3(WZ) < C 3 , 

inf po > p, and div_Bo = 0, then there is a solution (p,u,B) to (|l.l|) - (|1.5ll defined on R 3 x [0,T] 
satisfying 

p-p,B g C([0,T];iJ 3 (R 3 ))nC 1 ([0,T];iJ 2 (R 3 )) (1.28) 

and 

u E C([0, T];H 3 (R 3 )) n C 1 ([0, T};H 1 (R 3 )) n L 2 ([0, T]; iJ 4 (R 3 )). (1.29) 

T/ie equations in (|1.1[) - (|1.4[) are satisfied in the sense of equality of weak derivatives on R 3 x (0, T), each 
of these weak derivatives realized as an element ofC([0,T];H 1 (M 3 )), and the weak forms (|1.20j) - (|1.22j) 
hold. 

Furthermore, there is a positive number a depending on p, v, £, and P such that, if the above 
hypotheses hold with C3 < a, then the solution exists on all o/R 3 x [0, 00). 

As a by-product of our analysis we show that when the initial data (pq,uq, Bo) satisfies the 
smallness assumption (|1.23|) and is in H 3 , but is not necessarily small in H 3 , then the corresponding 
smooth solution in fact exists for all time. Since the weak solutions of Theorem 1.1 are constructed as 
limits of such global smooth solutions, we conclude that small-energy smooth solutions are dense in 
the set of small-energy weak solutions. The result is as follows: 

Theorem 1.2 Assume that the system parameters P,p,\ in (|1.1[) - (|1.4D satisfy the conditions in 
(1.13V ([1.16P and let 8 be as defined in (|1.15[) . Let positive numbers N and d < 5 be given. Then 
given initial data (pq,uq, Bq) satisfying (|1.17|) - (|1.19|l . the smallness condition (|1.23|) . and the regu- 
larity condition (po — P,uq,Bq) g H 3 (M. 3 ), the corresponding smooth solution (p,u,B) described in 
Kawashima 's theorem above exists on all o/R 3 x [0, 00). 
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This paper is organized as follows. We begin the proofs of Theorem 1.1 in section 2 with a number 
of a priori bounds for local-in-time smooth solutions. In section 3 we derive the necessary bounds for 
density by applying the estimates in Theorem 2.1 in a maximum principle argument along particle 
trajectories of the velocity, making important use of the monotonicity of P as described in (1.13). 
The small-energy assumption (|1 .23|) then enables us to close these arguments to show in Theorem 3.1 
that both the pointwise bounds for density and the a priori energy bounds of Theorem 3.1 do hold as 
long as the smooth solution exists. Finally section 4 we apply these now uncontingent a priori bounds 
to show that the H 3 norms of low-energy smooth solutions cannot blow up in finite time, and as a 
consequence that such solutions can be extended to all of {t > 0}, thereby proving Theorem 1.1 and 
1.2. 

We make use of the following standard facts (see Ziemer [24] Theorem 2.1.4, Remark 2.4.3, and 
Theorem 2.4.4, for example). First, given r € [2, 6] there is a constant C(r) such that for w £ if 1 (R 3 ), 

\\w\\ L r m < C(r) (iMlgTogf r ||V^||i 3 2 r (R ?/ 2r ) • (1.30) 

Next, for any r € (1, oo) there is a constant C(r) such that for w € W rl ' r (R. 3 ), 

IMU°°(R3) < C{r) \\w\\ w i, r (^3) (1.31) 

and 

(w)^ 3 <C(r)\\Vw\\ L r m , (1.32) 

where a — 1 — 3/r. Finally if V is the fundamental solution of the Laplace operator on R 3 , then there 
is a constant C such that for any / G _L 2 (M 3 ) n L 4 (R 3 ), 

||r*, * /]U~(R3) < C (H/IU^s) + ll/IUW . (1.33) 
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2. Energy Estimates 



In this section we derive a priori bounds for smooth solutions (p — p,u,B) 6 C([0,T}; H 3 (M. 3 )) of 
(|l.ip - (|1.4[) as described in Kawashima's theorem in section 1. These bounds will depend only on the 
quantities Co and d appearing in (jl . 1 7|) - (II . 19[) and will be independent of the initial regularity and 
the time of existence. Specifically, we define a functional A(t) for a given such solution by 

A(t)= sup (\\{p-p,u,B)(;s)\\% 2 +\\(p uUu B t (; S )\\ 2 L2 )+ f (||V U (-, S )||^ + |M-, S )|| 2 H1 )^, 

0<s<t JO 

(2.1) 

and we obtain the following a priori bound for A(t) under the assumptions that the initial energy Co 
in (|1.19j) is small and that the density remains bounded above and below away from zero: 

Theorem 2.1 Assume that the system parameters in (|l.l|l - (|1.4|) satisfy the conditions in (1.13)- ljl.l6|) 

and let 8 be as defined in (|1.15|) . Let positive numbers N and d < 6 be given. Then there are positive 
constants a, M , and 9 depending on the parameters and assumptions in (1.13)- ljl.l6|) and on a positive 
lower bound for d, such that: if (p,u,B) is a solution of () 1 . 1 1) - (1 1 .4[) on M 3 x [0, T] in the sense of 
Kawashima's theorem with initial data (po — P, Uo, Bo) € H 3 (M. 3 )) satisfying (|1.17I) - (I1.19I) with Co < a, 
and if 

P < p{x,t) <p onR 3 x [0,T], 

then 

A(T) < MC B . 



The proof will be given in a sequence of lemmas. It will be seen that the assumed regularity 
(|1.28l) - (jl.29l) suffices to justify the estimates that follow. The most subtle part in the analysis is to 
extract the dissipative structure of the system, which can be partially accomplished by introducing 
auxiliary variable functions w and v in Lemma 2.4. The methods we apply here are inspired by those 
of Lei-Liu-Zhou [15] for incompressible viscoelastic fluids. 

We begin with the following L 2 energy estimate: 

Lemma 2.2 Assume that the hypotheses and notations of Theorem 2.1 are in force. Then 

T + l {G(p)+ l - P \u\ 2 + l -\B\ 2 )dx+ [ (p\\7u\ 2 + \(dWu) 2 )dx = 0. (2.2) 

Proof. We multiply the momentum equation by u J , sum over j and integrate to obtain that for 
< t < T, 



t 

2 



^p\u\ 2 dx 



dx+ [ [ u ■ div [(i|£| 2 )/ 3x3 - BB T ]dxds 
JoJs. 3 



+ [[ u-VPdxds+ [ [ [p\Vu\ 2 + - p)(divu) 2 ] dxds = Q, (2.3) 



where the divergence of a matrix is taken row-wise. Next we define 

G = p f P S - 2 (P(s) - P)ds 
Jp 

where P = P(p), so that by the mass equation, 

G{p) t + div(G(p)u) + (P(p) - P)divu = 0. 
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Integrating and adding the result to (|2.3p we then get 

/ [\p\u\ 2 + G{p)] dx\l + f f u ■ div [(±\B\ 2 )I 3x3 -BB T ]dxds 
Jr 3 JqJr 3 

+ 11 [/i|Vu| 2 + ( div u) 2 ]dxds = 0. 

Finally we multiply the magnetic field equation by B and integrate to get 



/ i|B| 2 dx|*+ / / B ■ div (Bu T - uB T ) dxds = 0. 

JR 3 JoJrs 



(2.4) 



(2.5) 



We then obtain (|2.2|) by adding (I2.5[) to (12.4j) and using the fact that 

[u ■ div[(i|B| 2 )/ 3X 3 - BB T ] + B ■ div(Bu T ~ uB T )} dxds = 0. 

OJR 3 



Next we derive preliminary bounds for u, B and p in L°°([0, T]; 7? 2 (M 3 )) 



□ 



Lemma 2.3 Assume that the hypotheses and notations of Theorem 2.1 are in force. There is a positive 
number 6\ depends only on p and A such that for t 6 (0, T], 



d_ 

dt 



I (|Vii| 2 + p\u\ 2 + \VB\ 2 + \AB\ 2 + \B t \ 2 + \\7p\ 2 + \Ap\ 2 + \p t \ 2 )dx] + f {p\u\ 2 + \Vv 
Jr 3 J Jr 3 

I (p\u\ 2 + \Vii\ 2 + \\7u\ 2 + \AB\ 2 + \Ap\ 2 )dx +At[ \Vu\ 2 dx 
Jr 3 Jr 3 



7u\ 2 )dx 



<M\\(p-p,u,B)(;t)\\ 



01 

H 2 



(2.6) 



Proof. We multiply the momentum equation by u J , sum over j and integrate to obtain 

/ p\u\ 2 dx= I [-uVP + pu ■ Au + Am • V(divu) — u ■ V(^|-B| 2 ) + u ■ div(BB T )]da 
Jr 3 Jr 3 2 

Integrating by parts and absorbing terms, 



d_ 

dt 



\Vu\ 2 dx 



+ / p\u\ 2 dx < M 



(\VP\ 2 + \Vu\ 2 + \Vu\ 3 + \B\ 4 )dx 



(2.7) 



Next, we take the material derivative on the momentum equation, multiply it by ii 3 and integrate 

-— I p\ii j \ 2 dx= I -v? [P jt + div (Pju)]dx 
2 dt J R3 J R3 



pv? 



Au[ + div(Au j u) 



dx 



\o~ 5 u J [div(u)jt + div(div(u)jit)] dx 



dt 



+ I u J [-{±\B\ 2 ) jt - div((i|B| 2 ) iM )] dx + I v? [Aw{BB 3 ) t Aw{dw{BB J )u)} dx. 
Jr 3 Jr 3 

it 

p\u\ 2 dx 



Integrating by parts and absorbing terms 
d 



I \Vii\ 2 dx<M I {p 2 \X7u\ 2 + \B t \ 2 \B\ 2 + \Vuf)di 
Jr 3 Ur 3 

I p 2 \\7u\ 2 dx+ I iVu^dx 
Jr 3 Jr 3 



< M 



M\\B(;t)\\ 



I {\\7u\ 2 \B\ 2 + \VB\ 2 \u\ 2 )dx 
Jr 3 



(2.8) 



where the last inequality follows from the induction equation in fj 1 . 1 [) - fjl .4[) . 



For the term — 
dt 



integrate to obtain 



(\p t \ 2 + \B t \ 2 dx) 

/ \ P t\ 2 dx < [ \(pu) t \\Vp\d. 
at Jr 3 Jr3 
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, we differentiate the mass equation with respect to t and 



< / (\div(pu)\\u\ + \p\\u\ + \p\\Vu\\u\)\Vp\dx 



(2.9) 



For the terms / (| V£> + 1 AB\ )dx and / (|Vp + \Ap\)dx, we differentiate the mass equation 

JR 3 i JR 3 

with respect to x and integrate 

~f (\\7p\ 2 + \Ap\ 2 )dx < [ \Vu\\Vp\ 2 dx + [ \Vu\\p\\Vp\dx (2.10) 

2 dt J R 3 J R 3 J R 3 

+ / \D 3 x u\\p\\Vp\dx+ [ \Au\\V P \\Ap\dx+ [ \\7u\\Ap\ 2 dx 
Jr 3 Jr 3 Jr 3 



0i 



<M\\(p-p,u,B)(.,t)\\H 



f (\\7u\ 2 + \Au\ 2 + \D 3 x u\ 2 + \Ap\ 2 )da 
Jr 3 



Similarly, we differentiate the induction equation and integrate to obtain 

ld_ 

2dt 



(\VB\ 2 + \AB\ 2 )dx < M\\(p-p,u,B)(;t)\\% 



f (\\7u\ 2 + \Au\ 2 + \D 3 x u\ 2 + \A P \ 2 )dx 
Jm 3 



Adding (I2^- (I2~TT1) . 

/ (\\7u\ 2 + p\u\ 2 + \VB\ 2 + \AB\ 2 + \Vp\ 2 + \Ap\ 2 )dx]+ [ (p\u\ 2 + \\7ii\ 2 )dx 
Jm. 3 J Jm 3 



(2.11) 



(2.12) 



<M\\(p-p,u,B)(-,t)\ 



H' 2 



[ {\Au\ 2 + \D 3 x u\ 2 + \\7u\ 2 + \AB\ 2 + \A P \ 2 )dx 

Jm 3 



M / \Vu\ 2 dx. 



It remains to estimate the terms / |Au| 2 (fa; and / \D 3 u\ 2 dx. Using the definition of F and the 

Jr 3 Jr 3 



momentum equation, 

(p + \)Au 3 = F Xj +(p + XPlf + P Xj , 
AF = div( 5 ), 

where gi = pu? + (^\B\ 2 ) X . — div(-B-'-B). Therefore, by standard elliptic theory we obtain 



||Au(.,t)||£ a < M [\\WF(;t)\\h + \\Vu(;t)\\h + l|VP(-,*)||£ a + \\\VB\\B\(.,t)\ 



L 2 } ' 



||l£u(-,*)||£ a < M [||AF(.,i)||| 2 + ||Aa;(.,t)||| 2 + ||AP(,i)||| 2 + ||V(|VS||B|)(-,t)||| 2 ] 



||VF(-,*)||| 2 < M [\\pu(;t)\\h + |||V5||B|(-,i)| 
Using f2~T5l) - ([2~T7l) on (f2~T2"l) . (j2~6l> follows. 



(2.13) 
(2.14) 

(2.15) 

(2.16) 

(2.17) 
□ 



We extract the dissipative structure of the system by introducing auxiliary variable functions w 
and v^' as follows: 
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Lemma 2.4 Assume that the hypotheses and notations of Theorem 2.1 are in force. There is a positive 
number 82 depends only on p and X such that for t £ (0, T], 



d_ 
~dt 



/ (M 2 + yV j) | 2 )da; + / (|Vw| 2 + V \Vv^\ 2 )da 

jr3 f-j y K 3 f-j 



<M||(p-p,«,B)(.,t)| 



|0 2 

/f2 



/ {p\u\ 2 + |Vu| 2 + |Vu| 2 + |AS| 2 + |Ap| 2 )da 

JR3 



(2.18) 



where w = p 2 u + Vp and v^> = p 2 ii + Vi? 3 /or j = 1, 2, 3. 

Proof. We give the proof for w, the argument for is similar. By the definition of w, 

w t = Aw - A(p 2 u + Vp) + (p 2 ii + Vp) t , 
we then multiply the above by w and integrate to obtain 

\w\ 2 dx 



ld_ 

2 dt 



/ \Vw\ 2 dx = - A(p 2 ii + Vp) -wdx + (p 2 u + Vp) t -w. (2.19) 

Jr 3 JR3 JR3 

The first term in the right side of (|2.19[) is bounded by 

1 1 
\Vw\\V(p 2 u + V P )\dx < M (^J IV^I 2 ^ 2 (J (|p| 4 |Vw| 2 + |p| 2 |Vp| 2 |w| 2 + |Ap| 2 )dx^ 2 , 

(2.20) 



and the second term on the right side of ()2.19|) is bounded by 



|p*IIHIH 



{pu)tpwdx 



\pt\\Vw\dx 



(2.21) 



</ \&\v{pu)\\pu\\w\+ |V(^)|(|P t | + |Vw t | + |B t ||B|)da 



\p t \ 2 dx 



\Vw\ 2 dx 



where the last inequality follows by the mass equation and the momentum equation. Using (|2.20[1 - 
(|2.21[) on (12.19[) and absorbing terms, 



d_ 

dt 



l dx 



\Vw\ 2 dx 



<M\\(p-p,u,B)(;t)\\% 



[ (p\u\ 2 + \\7ii\ 2 + \\7u\ 2 + \AB\ 2 + \Ap\ 2 )dx 

Jr 3 



□ 



Proof of Theorem 2.1: First by the definition of w and 

3 

{\AB\ 2 + \Ap\ 2 )dx< [ (\p\ 4 \Vu\ 2 + \p\ 2 \Vp\ 2 \u\ 2 )dx+ [ (\\7w\ 2 + f2\\7v^\ 2 )dx. 
1 Jr 3 Jr 3 - =1 
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Therefore we add 



d_ 

dt 



< 



, (|2.18[) together and use the assumption p < p < p to obtain 
/ (\p - p\ 2 + \u\ 2 + \B\ 2 + \Vp\ 2 + \\7u\ 2 + \\7B\ 2 + \u\ 2 + \AB\ 2 + \Ap\ 2 )dx 

JS. 3 

3 

+ / (|Vm| 2 + H 2 + |Vm| 2 + |Vw| 2 + V|Vv w | 2 )dx 
Jr 3 ~1 

r 3 
/ (|Vu| 2 + \u\ 2 + |Vu| 2 + |Vw| 2 + V \Vv U) \ 2 )d. 
Jm 3 ~1 



M\\(p-p,U,B)(;t)\\H* 



3 = 1 



3=1 

+ M \Vu\ 2 dx, 
(2.22) 



for some M > and #3 > 0. We define a functional X(t) by 

X(t)= sup / (\p- p\ 2 + \u\ 2 + \B\ 2 + \\7p\ 2 + |Vit| 2 + |VB| 2 + \ii\ 2 + \AB\ 2 + \Ap\ 2 )(x, s)dx 

0<s<t JR3 

+ ff (\Vu\ 2 + \ii\ 2 + \Vii\ 2 + \Vw\ 2 + \Vv {l) \ 2 ){x,s)dxds, 



oJr 3 



then 



SUP \\(p-p,U,B)(;s)\\ 2 H3 <X(t), 
0<s<t 



and so by integrating (|2 . 22|) from to t, we get 

X(t) < MX(t)^ +1 +X(0) + M [ [ \Vu\ 2 dx 

< MX{t)^ +1 + X(0) + Mp^Co, 

where the last inequality follows from (|2.2[) . Using the fact that X(t) is continuous in time, there are 
positive constants a,M, and 9 such that if Co < a, then 



A(T) < X(T) < MC$, 



which finishes the proof of Theorem 2.1. 



(2.23) 

□ 



3. Pointwise bounds for the density 

In this section we derive pointwise bounds for the density p, bounds which are independent both of time 
and of initial smoothness. This will then close the estimates of Theorem 2.1 to give an uncontingent 
estimate for the functional A defined in (12.11). The result is as follows: 



Theorem 3.1 Assume that the system parameters in (|l.ll) - (jl.4l) satisfy the conditions in (1.7)- (|1.16|) 

and let positive numbers d < S be given. Then there are positive constants a, M, and 8 depending 
on the parameters and assumptions in (1.7) - (|1.16|) and on a positive lower bound for d, such that: if 
(p,u,B) is a solution of (|l.ip - (|1.4j) on R 3 x [0,T] in the sense of Kawashima's theorem with initial 
data (p - p,u ,B ) 6 H 3 (R 3 )) satisfying (f!~ll^ - (fi~T9|) with C Q < a and p(x,t) > on R 3 x [0,T], 
then in fact 



and 



P < p{x,t) < p on R 3 x [0,T] 



A(T) < MC$. 



(3.1) 



(3.2) 
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The proof of Theorem 3.1 consists of a maximum-principle argument applied along integral curves 
of the velocity field. We begin with two auxiliary propositions: 

Lemma 3.2 Let (p,u,B) be as in Theorem 3.1 and suppose that < c-i < p < c\ onl 3 x [0, T]. Fix 
to > and define the particle trajectories x : [0, oo) x K 3 4 I 3 by 

( x(t,y) = u(x(t,y),t) 
\ x{t ,y) = y. 

Then there is a constant C depending only on c\ and C2 such that if g G L 1 (M 3 ) is nonnegative and 
t G [0,T], then each of the integrals / g(x(t,y))dy and / g{x)dx is bounded by C times the other. 

Proof. The proof is exactly as in Hoff [6], pg. 27. □ 



Proof of Theorem 3.1: First we choose positive numbers b and b' satisfying 

p<b<p + d<p-d<b'<p. 

Recall that po takes values in [p + d, p — d] , so that p G [p, p] on R 3 x [0, r] for some positive r by the 
time regularity (|1.28j) . It then follows from Theorem 2.1 that A(t) < MC®, where M is now fixed. 
We shall show that if Co is further restricted, then in fact b < p < b 1 on R 3 x [0,r], and therefore by 
a simple open-closed argument that b < p < b' on all of R 3 x [0, T], and therefore that A(T) < MCq 
as well. We shall prove the required lower bound, the proof of the upper bound being similar. 
Fix y G R 3 and define the corresponding particle path x(t) by 

x{t) = u(x(t),t) 
x(0) = y. 

Suppose that there is a time t\ < r such that p{x{t\), t\) = b. We may take t\ minimal and then 
choose to < ti maximal such that p(x(to),to) = p + d. Thus p(x(t),t) G [b,p + d] for t G [io,ti]. We 
consider two cases: 

Case 1: t < t t < T A 1. 

We have from the definition (jl,12|) of F and the mass equation that 

Pj t [log p(x(t),t) - log(p)] + P(p(x(t),t)) -P= -F(x(t),t). 
Integrating from to to t\ and and abbreviating p(x(t),t) by p(t), etc., we then obtain 

p[log p(s) - log(p)] U + f 1 [P(s) - P]ds = - f 1 F(s)ds. (3.3) 

to Jt Jt 

We shall show that 

F(s)ds < MC 9 o (3.4) 

to 

for a constant M which depends on the same quantities as the M from Theorem 2.1 (which has been 
fixed). If so, then from (|3.3jl . 

//[log b - log(p + d)]>- / [P(s) - P]ds - MC e o > -MC 9 o, (3.5) 



where the last inequality holds because p(t) takes values in [b, (p + d)] C [p, min{p, p'}] for t G [to, ti], 
and P is increasing on [p, p'] (see (1.7)- (|1.15|) 'l. But ()3.5|) cannot hold if Co is small depending on 
M , b, and p + d. Stipulating this smallness condition, we therefore conclude that there is no time t\ 
such that p(t\) = p(x(ti),ti) — b. Since y G R 3 was arbitrary, it follows that b < p on R 3 x [0, r], as 
claimed. The proof that p < b' is similar. 
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To prove (|3.4p we recall the definition of F and apply Theorem 2.1 to obtain 



F(s)ds< / \\F(;S 



: ds 



< 



< 



(jM + \)\\dxVu(;a)\\ L oo + \\(P-P)(.,8)\\ 



ds 



(M + A)(||Vu(- ) »)|| i * + ||Att(.,*)|| i + ll(i'--P)(-.«)IU* + l|VP(-,«)|| £ 4 
< MA(T) < MC e . 



ds 



Case 2: 1 < < < *i- 

Again by the mass equation and the definition (|1.12[) of F, 

j t { P (t) -p)+ fi-'pmpit) - p(p)) = -p- i P {t)F{t). 

Multiplying by {pit) — p) we get 

\j t {p{t) P? + /x-VWpWCpW - P? = -p-'pitmt) - p)F{t) (3.6) 



where /(*) = (P{t) - P{p)){p{t) - p) -1 . Since f(t) > f(t ) > on [t ,*i], we have that 

M _1 /(*)P(*)(P(*) - P) 2 > »- X f(to)p(p(t) - pf for i e [to.ti] 
and therefore from (13.61) that 



(p( tl ) - p) 2 - (p(t ) - p) 2 < M r \\F{;s)\\ 2 Laa ds. (3.7) 

Jto 



We shall show that 

2 



||F(-,a)||i.»cfa <MC B , (3.8) 



to 



so that from (|3.7|) . 

< (6 - p) 2 - (p + d - p) 2 < MC e Q . 

This cannot hold if Co is sufficiently small, however, so that, as in Case 1, there is no time t\ such that 
p(t\) = p(x(ti),ti) = b. Since t/6R 3 was arbitrary, it follows that b < p on R 3 x [0, r], as claimed. 
To prove ([III) we a PP!y <P~^2l> and ([L33} to get 

f 1 \\F(;s)\\ 2 Lao ds< f 1 [\\g(;s)\\ 2 Li + \\g(;s)\\ 2 L2 ]ds. (3.9) 

The first integral on the right side above is bounded as follows: 



! 1 ||ff(-,s)|| i4 d S < / 1 [\\pu{;s)\\ 2 Li + \\WB.B(.,s)\\ 2 Li ] ds 

Jt Jt 

rti < 

< / [H-, S )||| 2 ||Vu(-, s)\\h + IIVB • B(-, s)\\l 2 \\D x (VB ■ B)(;s)\\\ 

Jt L 

t ~ t ~ 

[ 1 [ \u\ 2 dxds] ( [ 1 [ \Vu\ 2 dxds ) 
A Jr 3 / \Ji JR 3 / 

/ / \VB ■ B\ 2 dxds ) ( / / \D X (VB ■ B)\ 2 dxds j <MC O: 
Ji Jr 3 / \Ji Jr 3 / 



where the last inequality follows from Theorem 2.1. The second integral on the right side of (J3T9J) can 
be bounded in a similar way, and Q3.8P is proved. □ 
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4. Global Existence of Smooth Solutions: Proof of Theorem 1.1 and 1.2 

In this section we prove the global-in-time existence of smooth solutions with initial data satisfying 
the condition (|1.23[) . Specifically we show in Theorem 4.1 below that, as a consequence of the time- 
independent bounds p.ip and (I3.2p in Theorem 3.1, the ff 3 (R 3 ) norms of such solutions remain finite 
in finite time. Theorem 1.1 and 1.2 then follow immediately. 

Theorem 4.1 Assume that the system parameters in (| 1 . 1 [) - fjl .4[) satisfy the assumptions made in The- 
orem 1.1 and let d, 5, a, 9, and M be as described in Theorem 3.1. Then there is a continuous function 
M'(T) such that, if (p,u, B) is a smooth solution of (|l.l[) - (|1.4p on R 3 x [0, T] as described in Theorem 
3.1 with initial data (po — P, uq, Bq) S H 3 (M. 3 ) satisfying (|1.17l) - (|1.19[) as well as the condition (|1.23[) . 
then 

sup \\(p-p,u,B)(;s)\\ H3 + [ \\u(; S )\\ 2 H4 ds<M'(T). (4.1) 



o<s<t Ja 

Proof. We give the proof in a sequence of steps which are similar to those given in Suen and Hoff [23 
pp. 26-31. Throughout this proof we let M' denote a generic function of T as described above. 

Step 1: We record the following bounds for spatial derivatives of u in terms of its time derivative and 
lower-order terms: 

\\D 3 x u(;t)\\ L 2 < M [||Vp.u(-,t)|| ia + \\p\7u(;t)\\ L 2 + \\B- DlB(;t)\\ L * 

+ |||V J B| 2 (.,t)|| L2 + ||^P( /J )(.,t)|| L2 ]. (4.2) 

These follow immediately from the momentum equation and the ellipticity of the Lame operator 
L = pA + AVdiv. □ 

Step 2: The velocity vector satisfies the following bounds: 

sup ||Vu(-,i)IU" + / / \D 2 x u\ 2 dxdt<M'{T) (4.3) 

0<t<T Jo JR3 

and 

sup \\u(;s)\\ H s <M'(T). (4.4) 

0<s<T 

Proof: First for h > we define the forward difference quotient D' t l by 

*c/)w- /( ' + *>- /w . 

and we let $ = D^{u^) + u ■ Vu-? . Then by differentiating and differencing the momentum equation 
and applying the bounds in (13. 1^ - (|3.2|) in an elementary way, we obtain 

o\(f )xj \ 2 dx+ [[ (|V^ 3 .| 2 + |D t h (divK.)) +u- V{div(u Xj ))\ 2 )dxds 

<M'(T)+ I I \Vcj)\ 2 dxds + 0(h), 

where O(h) — > as h — > and ()4.3|) follows immediately. We then substitute the results in (|3.ip - (|3.2|) 
and (1431) into (l4~2l) to obtain the H 3 bounds in (Oil. □ 



Step 3: The following preliminary bounds hold for the L 2 ([0, T); H 4 (M. 3 )) norms of u: 

\D x u\ 4 dxds < M'(t) + M I I (\D 3 x p\ 2 + \D 3 X B\ 2 ) dxds. (4.5) 



Proof: These bounds follow immediately by differentiating the momentum equation twice with respect 
to space, expressing the above fourth derivatives of u in terms of second derivatives of ii, Vp, VB, and 
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lower order terms via the ellipticity alluded to in Step 2, and then applying the bounds in 
and (HU). □ 



Step 4: The following bounds for p, u and B hold: 

sup (\\D 3 x p(;t)\\L* + \\DlB(;t)\\ L 2)+ [ [ \D 4 x u\ 4 dxds < M' (T) . (4.6) 

0<t<T JoJR 3 

Proof: We apply two space derivatives and one spatial difference operator Dj to the mass equation 
to obtain an evolution equation for D^D Xi D Xk p, then multiply by this same quantity and integrate. 
Proceeding exactly as in Step 4 and applying (|4.5j) we obtain 

/ \D^D Xt D Xk p\ 2 dx<M'(T) + f f \D*D Xt D Xk p\ 2 dxds + M f f \D 4 u\ 2 dxds 

< M'(T) + M [ { \D 3 x p\ 2 dxds. 
JoJr 3 

Wc let h — > and apply Gronwall's inequality to obtain the required bound for the first term on 
the left in (14.61) . The bound for B can be obtained in a similar way. Then substituting into (|4.5[) 
we obtain a bound for the second. This proves (|4.6|) and together with (14. 3} completes the proof of 
Theorem 4.1. □ 

Proof of Theorem 1.1 and 1.2: It suffices to prove Theorem 1.2 since Theorem 1.1 follows immediately 
by standard compactness argument. The proof of Theorem 1.2 involves an open-closed argument on 
the time interval which is identical to the one given in Suen and Hoff 23J pp. 31. We omit the details 
here. □ 
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